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Ferrogels and magnetic elastomers differentiate themselves from other materials by their unique 
capability of reversibly changing shape and mechanical properties under the influence of an exter¬ 
nal magnetic field. A crucial issue in the study of these outstanding materials is the interaction 
between the mesoscopic magnetic particles and the polymer matrix in which they are embedded. 

Here we analyze interactions between two such particles connected by a polymer chain, a situation 
representative for particle-crosslinked magnetic gels. To make a first step towards a scale-bridging 
description of the materials, effective pair potentials for mesoscopic configurational changes are spec¬ 
ified using microscopic input obtained from simulations. Furthermore, the impact of the presence of 
magnetic interactions on the probability distributions and thermodynamic quantities of the system 
is considered. The resulting mesoscopic model pair potentials can be used to economically model the 
system on the particle length scales. This first coarse-graining step is important to realize simplified 
but realistic scale-bridging models for these promising materials. 


I. INTRODUCTION 

Ferrogels and magnetic elastomers are fascinating ma¬ 
terials, born by the union of polymeric networks and fer- 
rofluids. Their amazing properties derive from the unique 
combination of the elastic behavior typical for polymers 
and rubbers [1] on the one hand, and magnetic effects 
characteristic of ferrofluids and magnetorheological fluids 
[2-10] on the other. One of the most interesting results is 
that their shape and mechanical properties can be exter¬ 
nally controlled by applying a magnetic field [11-21]. A 
form of tunability, distinguished by reversibility as well 
as non-invasiveness and based on a magneto-mechanical 
coupling is one of the most appealing properties of these 
materials. This makes them excellent candidates for 
the use as soft actuators [22], magnetic field detectors 
[23, 24], as well as tunable vibration and shock absorbers 
[14, 25]. Moreover, studying their heat dissipation due 
to hysteretic remagnetization in an alternating external 
magnetic field might be helpful to understand better the 
processes during possible applications in hyperthermal 
cancer treatment [26, 27]. 

Typically, these materials consist of cross-linked poly¬ 
mer networks in which magnetic particles of nano- or 
micrometer size are dispersed [16]. A central role in 
the coupling of magnetic and mechanical properties is 
played by the specific interactions between the embed¬ 
ded mesoscopic magnetic particles and the flexible poly¬ 
mer chains filling the space between them. These cou¬ 
plings are responsible for a modified macroscopic elas¬ 
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ticity [12, 28, 29], orientational memory effects [30, 31], 
and reversibility of the magnetically induced deforma¬ 
tions [15]. 

Many theoretical and computational studies have been 
performed on the topic, using different approaches to 
incorporate the particle-polymer interaction. Some of 
them rely on a continuum-mechanical description of 
both the polymeric matrix and the magnetic component 
[23, 32, 33]. Others explicitly take into account the dis¬ 
crete embedded magnetic particles, but treat the polymer 
matrix as an elastic background continuum [34, 35]. Usu¬ 
ally in these studies, an affine deformation of the whole 
sample is assumed. The limitations of such an approach 
for the characterization of real materials have recently 
been pointed out [36]. In order to include irregular distri¬ 
butions of particles and non-affine sample deformations, 
other works employ, for instance, finite-element methods 
[19, 37-40]. 

To optimize economical efficiency, a first step is the use 
of simplified dipole-spring models. In this case, steric re¬ 
pulsion and other effects like orientational memory terms 
can be included [30, 36, 41-44]. A step beyond the of¬ 
ten used harmonic spring potentials can be found in [42] 
where non-linear springs of finite extensibility are consid¬ 
ered. 

From all the studies mentioned above it becomes clear 
that microscopic approaches that explicitly resolve poly¬ 
mer chains are rare [45, 46], and what is particularly 
missing is a link between such microscopic approaches 
and the mesoscopic models that only resolve the mag¬ 
netic particles, not the single polymer chains. In partic¬ 
ular, a microscopic foundation of the phenomenological 
mesoscopic expressions for the model energies should be 
built up. 

The present work is a first step to close this gap. 
We consider an explicit microscopic description in a 
first simplified approach: a single polymer chain, dis¬ 
cretized through multiple beads, each representing a 
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coarse-grained small part of the polymer, connects two 
mesoscopic particles. The ends of the polymer chain 
are rigidly anchored on the surfaces of the two meso¬ 
scopic particles, which are spherical, can be magnetized, 
and are free to rotate and change their distance. From 
molecular dynamics simulations on the microscopic level, 
we collect the statistics of the micro-states correspond¬ 
ing to different configurations in the mesoscopic model. 
Based on these statistics we derive effective mesoscopic 
pair potentials, inspired by previous achievements for 
other polymeric systems [47, 48]. The subsequent step 
in scale-bridging, connecting the mesoscopic picture to 
the macroscopic description, has been recently addressed 
[49] for a special class of magnetic polymeric materials. 

In the following, we first define and describe our model 
in section II. Then, in section III, we mention the de¬ 
tails of the microscopic simulation. In section IV, we 
further characterize the probability distribution of the 
mesoscopic variables, connecting it to a wrapping effect 
in section V. After that, in section VI, we determine 
the values of mesoscopic model parameters based on the 
results of our microscopic simulations. In section VII, 
we derive an approximated expression for the mesoscopic 
effective pair potential characterizing the particle con¬ 
figurations. In this way, we build the bridge from the 
explicit microscopic characterization to the mesoscopic 
particle-resolved models by averaging over the micro¬ 
scopic details. Last, in section VIII, we consider the ef¬ 
fect of adding magnetic moments to the particles and, 
in section IX, the effect of increasing magnetic interac¬ 
tion on the thermodynamic properties, before we draw 
our final conclusions in section X. Appendix A addresses 
the trends in the dependences of the mesoscopic model 
parameters on varying microscopic system parameters, 
while appendix B briefly comments on the separability 
of magnetic interactions in the mesoscopic picture and 
microscopic chain configurations. 


II. THE SYSTEM 


Our simplified system is composed of two mesoscopic 
and spherical particles, both of radius a, and a polymer 
chain explicitly resolved by N = 60 beads of diameter a 
and interconnected by harmonic springs. Here we choose 
the mesoscopic particle radius a to be 5 a. 

We consider steric repulsion between all described 
particles through a WCA potential, which represents a 
purely repulsive interaction. It is given by 


Vwca 




for r' < r c 


0 otherwise 


(i) 

where r' is the distance between the particle centers, 5 
denotes the energy scale of the potential, and r c = 2 1 / 6 cr / 
is the cut-off distance. For any combination of large and 


small particles, a' is chosen as the sum of the radii of 
the respective particles, which is equivalent to the mean 
of their respective diameters. In our simulations, we set 
e = IO&bT, where is the Boltzmann constant and T is 
the temperature of the system. Neighboring beads within 
the chain are bound by means of a harmonic potential 

V H (r') = I k(r' - r 0 ) 2 , (2) 

where we choose the force constant k = 10 ksT/a 2 . The 
equilibrium distance ro is set to match the cut-off dis¬ 
tance of the WCA-potential ro = r c = 2 1 / 6 cr. 

The ends of the chains are bound via the same har¬ 
monic potential to binding sites placed below the surface 
of the mesoscopic particles, see figure 1. These bind¬ 
ing sites are rigidly connected to the mesoscopic parti¬ 
cles and follow both, their translational and rotational 
motion. Thus, when the magnetic particle moves or ro¬ 
tates, the binding site of the polymer chain has to follow, 
and vice versa. The technical details for the virtual sites 
mechanism can be found in [50]. We identify the an¬ 
choring points of the polymer chain by the vectors a\ 
and a 2 , respectively (see figure 1). Furthermore, the dis- 



Figure 1: A simplified sketch of the geometry of the micro¬ 
scopic system, r is the vector connecting the centers of the 
mesoscopic particles, ai and a 2 identify the anchoring points 
of the polymer chain on the surfaces of the particles. 6 1 , O 2 , 
and (j) are the angles that represent the remaining relative 
rotational degrees of freedom of the system. 

tance vector between the two mesoscopic particle centers 
is indicated as r, with magnitude r. The angles that the 
vectors oq and a 2 adopt with respect to the connecting 
vectors r and -r are denoted as 6 \ and 62 , respectively. 
01 and 62 are zenithal angles and they can span the in¬ 
terval [0,7r]. Last, cj) G [—7r,7r[ is the relative azimuthal 
angle between the projections of cq and (12 on a plane 
perpendicular to r. It can be used to parametrize the rel¬ 
ative torsion between the two particles around the r axis. 
Let us for brevity introduce the vector 7 = (r, 0i, # 2 , 0). 
Therefore, the 7 -space, where our mesoscopic vector 7 is 
defined, is given by [0, +oo[x [0, tt\ x [0,7r] x [—7r,7r[. 
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Through molecular dynamics simulations (thoroughly 
described in the next section) we find the probability 
density p c (l) of a certain configuration 7 of the two 
mesoscopic particles. We normalize such that 

f = 1 , with d 7 = sin( 0 i) sin(# 2 )drd#id# 2 d</> the 

7 -space volume element. Moreover, we are working in 
the canonical ensemble. 

From statistical mechanics [51], we know that the 
probability density p c (l) to find the system in a cer¬ 
tain configuration 7 is p c ( 7 ) = exp[—/3V C (7)]/Z c , where 
/3 = 1 /k B T and Z c is the partition function of the system. 
Shifting W( 7 ) by an appropriate constant, we can still 
reproduce p c {l) but simultaneously normalize Z c = 1. 
Then, since ln(Z c ) = 0, 

Vc{ 7 ) = -k B T In [p c ( 7 )] (3) 

represents an effective energy of the state 7 of our meso¬ 
scopic two-particle system and corresponds to an effec¬ 
tive pair potential on the mesoscopic level, see also [47]. 
Through the normalization of Z c , we set our reference 
free energy F = — k B Tln(Z c ) equal to zero. 

III. MICROSCOPIC SIMULATION 

To obtain the probability distribution from which the 
mesoscopic pair potential is derived, we performed molec¬ 
ular dynamics simulations using the ESPResSo software 
[50, 52]. Because entropic effects are important to cap¬ 
ture the behavior of the polymer, the canonical ensemble 
is employed. This is achieved using a Langevin ther¬ 
mostat, which adds random kicks as well as a velocity 
dependent friction force to the particles. For the trans¬ 
lational degrees of freedom of each particle, the equation 
of motion is then given by 

m p v(t) = ~(v(t) + F r + F, (4) 

where m v is the mass of a particle, F is the force due 
to the interaction with other particles, F r denotes the 
random thermal noise, and ( is the friction coefficient. 

To maintain a temperature of T, according to the 
fluctuation-dissipation theorem, each random force com¬ 
ponent must have zero mean and variance 2 k B T(. Fur¬ 
thermore, random forces at different times are uncorre¬ 
lated. In order to track the orientation of the magnetic 
nanoparticles, rotational degrees of freedom also have to 
be taken into account. This is achieved by means of a 
Langevin equation of motion similar to (4) where, how¬ 
ever, mass, velocity, and force are replaced by inertial 
moment, angular velocity, and torque, respectively [53]. 

The time step for the integration using the Velocity- 
Verlet method [54] is dt = 0.01. To sample the proba¬ 
bility distribution, we record the state of the simulation 
7 every ten time steps. In order to obtain a smooth 
probability distribution over a wide range of parameters, 


34 billion states have been sampled in total, by running 
many parallel instances of the simulation, summing up 
to about 10 4 core hours of CPU time. 

Finally, we find the probability distribution by sort¬ 
ing the results of our simulations into a histogram with 
100 bins for each variable (r, 6 1 , # 2 , and 0). When a 
64-bit unsigned integer is used as data type, this leads 
to a memory footprint of 800 Mb. Hence, the complete 
histogram can be held in memory on a current computer. 
If the resulting numerical version of the effective pair po¬ 
tential as defined in (3) is to be used in a simulation, a 
smoothing procedure should be employed, especially in 
parts of the configuration space with a very low proba¬ 
bility density. One approach here might be hierarchical 
basis sets. 

In our simulations, we chose the thermal energy k B T 
as well as the mass m p , the friction coefficient £, and the 
corresponding rotational quantities to be unity. We mea¬ 
sure all lengths in units of a and the energies in multiples 
of k B T. 


IV. DESCRIPTION OF THE PROBABILITY 
DENSITY 

We now consider some aspects of the probability den¬ 
sity resulting from the microscopic simulation described 
in section III. The entropic role of the microscopic de¬ 
grees of freedom is considered by assigning to every con¬ 
figuration 7 a certain probability, given by the number 
of times it was encountered in the simulation divided by 
the total number of recorded states. As explained be¬ 
fore, 7 = (r, # 1 , # 2 , (j)) corresponds to the set of variables 
that we use to describe the state of the system on the 
mesoscopic level. 

In calculating the probability density from the micro¬ 
scopic simulations, we must remember the normalizing 
condition 

p c (l) sin( 6 > i) sin( 6 > 2 )drd 6 >id 6 > 2 d 0 = 1, (5) 

where r is integrated over [0,+oo[, 0 1 and 62 over [0,7r], 
and <j) over [— 7 r, 7 r[. Therefore, to obtain the probability 
density, we have to properly divide the data acquired 
from the simulations by the factor sin( 6 > i) sin^). 

It is useful to introduce here the average of a quan¬ 
tity over p c , defined as (*) c = f • p c ( 7 )d 7 . We can 
therefore calculate the average value of 7 , 7 = (j) c , 
and the covariance matrix = ((a - a) (0 - /?)) , for 
a,/3 = r, # 1 , # 2 ,</>. We find 7 ~ (20.20cr, 0 . 367 T, 0 . 367 T, 0). 
It is more practical to discuss the system in terms of 
correlation than in terms of covar iance. Correlation is 
defined as Q c a p = pp (no summation rule 

in this expression), is dimensionless, and g c a p G [— 1 , 1 ]. 
Here, we obtain 
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Q 


c 


( 1 
-0.341 ±0.028 

-0.356 ±0.029 

\ 0.083 ± 0.022 


-0.341 ± 0.028 -0.356 ± 0.029 0.083 ± 0.022 ' 
1 -0.006 ±0.006 -0.083 ±0.010 

-0.006 ±0.006 1 -0.083 ±0.011 

-0.083 ±0.010 -0.083 ±0.011 1 


( 6 ) 


where lines and columns refer to r, #i, 62 , 0 i n this order. 
The diagonal elements are unity by definition, because 
each variable is perfectly correlated with itself. The er¬ 
rors follow from the unavoidable discretization during the 
statistical sampling procedure in the simulations, where 
the results have to be recorded in discretized histograms 
of finite bin size. 

We find a strong anticorrelation between r and (9 i ? 2, 
meaning that when the distance between the mesoscopic 
particles changes they tend to rotate. We will address 
in detail the background of this behavior in section V in 
the form of the wrapping of the polymer chain around the 
mesoscopic particles. For angles 6 \ and O 2 different from 
0 and 7 r it is clear that such a wrapping and correspond¬ 
ing distance changes are likewise induced by modifying 
the relative torsion of angle 0. These coupling effects are 
partly reflected by the remaining non-vanishing correla¬ 
tions which, however, are weaker than the correlations 
between r and # 1 , 62 by at least a factor 4. The cor¬ 
relations between 6 \ and 62 are very weak and, within 
the statistical errors, may in fact even vanish. Within 
the statistical errors, the magnitudes of the correlations 
Q c r 0i and g c r 02 as well as Q c 0i (j) and q % 2 ^ agree well with 
each other, respectively, which reflects the symmetry of 
the system. Finally, we performed additional microscopic 
simulations for different sizes of mesoscopic particles and 
different chain lengths of the connecting polymer. As a 
general trend, we found that the correlations tend to de¬ 
crease in magnitude for smaller mesoscopic particles and 
for longer polymer chains (see appendix A). 

As a further step in the analysis of p c , we can de¬ 
termine the marginal probability density p a (<a) for one 
of the four mesoscopic parameters a = r,0 1 , #2 , 0 in¬ 
tegrating out the other three, for instance, po 1 (0 1 ) = 
fPc{ 7 ) sin(# 2 )drd 02 d</>. This is the total probability den¬ 
sity for the variable 6 \ to assume a certain value, regard¬ 
less of the others. Calculations for p r (r), pe 2 {^ 2 ), and 
P 0 ( 0 ) are analogously performed by integrating out all re¬ 
spective other variables (see figure 2 ). Of course, p a (<a) is 
still a normalized probability density since f p a (a)da = 
f Pc (7)^7 = 1, where we indicate da = sin(^)d 6 ^ for 
a = 0i (i = 1,2) and da = dr, d0 for a = r, 0. Follow¬ 
ing (3), we also introduce the one-variable effective pair 
potentials 

V a (a) = -k B T\n\p a (a)} (7) 

that are associated with the corresponding single-variable 
marginal probability density. 

In figure 2 most of the probability density p r (r ) for 
the interparticle distance is contained between r = 10 a 



Figure 2: Marginal probability densities for the single vari¬ 
ables. The abscissa has been rescaled: r ranges from 0 to 
lOOcr, whereas the angles range from 0 to 7r. pe 1 (0i) and 
pe 2 {0 2 ) are practically indistinguishable due to the symmetry 
of the set-up. Due to the symmetry of the system under the 
transformation 0 —— 0, here we only plot 0^,(0) from 0 to 
7 r. The maxima of the single-variable densities are located 
at r — 17cr and 0 1,2 = 0 = 0, respectively. The maxima of 
pe 1 (# 1 ) and pe 2 (O 2 ) at # 1,2 — 0 are found from the microscopic 
simulations after the 7 -space normalization contained in (5) 
has been taken into account. 

and r = 50<j, with a single maximum at r = 17cr. More¬ 
over, the steep increase in p r (r) at r = 11 a is to be 
attributed to the WCA steric repulsion, since at that 
distance the two particles are in contact. From figure 2 
we find that the maximum of pe^Oi) (with i — 1 , 2 ) is 
located at 0i = 0. The highest probability density for 
#1,2 = 0 is obtained from the microscopic simulations by 
taking into account the 7 -space normalization following 
from the use of spherical coordinates in (5), see also [47]. 
Moreover, p^(0) shows a maximum for 0 = 0, indicating 
that, as expected, the system does not tend to sponta¬ 
neously twist around the connecting axis in the absence 
of further interactions. The presence of a maximum at 
0 = 0 confirms that p^ is an even function of 0 invariant 
under the transformation 0 —>> — 0 , as expected from the 
symmetry of the set-up. 

V. WRAPPING EFFECT 

Before developing an approximate analytical expres¬ 
sion for the effective pair potential between the meso¬ 
scopic particles, it is helpful to examine in detail the re- 
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suits of the molecular dynamics simulations. In a mag¬ 
netic gel in which mesoscopic magnetic particles act as 
cross-linkers [31, 55], two driving mechanisms for a defor¬ 
mation in an external magnetic field are possible. First, 
in any magnetic gel, the magnetic interactions between 
the mesoscopic particles lead to attractions and repul¬ 
sions between them, which directly implies deformations 
of the intermediate polymer chains. As we would like 
to examine this mechanism separately, the magnetic in¬ 
teraction was not included explicitly in the simulations 
described in section III. Rather, it will be considered 
later in section VIII. Second, due to the anchoring of the 
polymer chains on the surfaces, rotations of the meso¬ 
scopic particles are transmitted to chain deformations. 

It has been shown in model II of [45] that the second 
mechanism on its own can lead to a deformation of such 
a gel in an external magnetic field: if mesoscopic mag¬ 
netic particles are forced to rotate to align with the field, 
the polymers attached to their surfaces have to follow. 
The resulting wrapping of the polymer chains around 
the particles leads to a shrinking of the gel. Transfer¬ 
ring this to our model system, it could imply an external 
magnetic field rotating the particles to a state in which 
the angles 6 \ and O 2 are non-zero. The effect of induced 
particle rotations on the interparticle distance is illus¬ 
trated in figure 3, where the effective pair potential V c (r) 
is plotted for the case of </> = 0 and various values of 
61 = 62 = 0 . I.e., both particles are rotated by the 

same angle 6 with respect to the connecting vector =br, 
respectively. It can be seen that the more the particles 
are rotated, the more the minimum of the effective pair 
potential is shifted towards closer interparticle distances 
corresponding to smaller separation distances r. 

To get a more detailed picture, we also consider in¬ 
dependent rotations of the two mesoscopic particles. In 
Fig. 4, the average distance between the particles is de¬ 
picted in a contour plot as a function of the angles 0 \ 
and O 2 • Images are shown for torsion angles of 0 = 0 and 
4> = 7r. It can be seen that quite a decrease in the av¬ 
erage distance can already be achieved by rotating only 
one particle. However, very strong reduction in inter¬ 
particle distance can only occur when both particles are 
rotated. The torsion angle <p does not alter the general 
trend of reduction of the average distance when the par¬ 
ticles are rotated. However, the resulting numbers vary 
to a certain degree. 


VI. ONE-VARIABLE EFFECTIVE PAIR 
POTENTIALS 

We now introduce some mesoscopic analytical expres¬ 
sions to model the effective pair potentials V a (a) intro¬ 
duced in (7). We will determine functional forms and 
parameters that can be used to model strain and torsion 
energies. Here, we use the term ” strain“ to denote varia¬ 
tions of r, whereas the term ” torsion “ is used to describe 
changes in the angles 6 \ and 62 or 0, depending on the 



r/cr 


Figure 3: Plot of the effective mesoscopic pair potential V c (r) 
for a situation in which both particles are rotated by the same 
amount in the angles 6 \ — 62 = 0 at a torsional angle of 
<f> — 0. It can be seen that the more the particles are rotated 
out of their equilibrium position, the closer they will approach 
each other because the minimum of the effective pair potential 
shifts to smaller separation distances r. The irregularities in 
the effective pair potential for high values of r are attributed 
to the naturally low sampling of those low-probability config¬ 
urations. 


initial orientation of the spheres and relative rotations. 
The quality of the analytical model expressions will be 
analyzed by fitting to the corresponding results from the 
microscopic simulations. 

First, we turn to the energetic contributions arising 
from changes in the interparticle distance r. The corre¬ 
sponding effective energy V r (r) obtained from the micro¬ 
scopic data is plotted in figure 5. It can be seen that there 
are essentially two regimes: up to r ~ 11a the WCA re¬ 
pulsion between the two mesoscopic particles dominates, 
whereas, for r > 11 a, V r (r) shows a smooth behavior and 
a single minimum arising from the entropic contribution 
of the polymer chain. Moreover, at r > 55a, V r (r) shows 
an irregular, non-smooth behavior. This is attributed to 
the low sampling rate of this extremely stretched config¬ 
uration, which has a very low probability to occur in the 
microscopic simulations (see figure 2). As a first approx¬ 
imation, it is natural to reproduce V r (r) by a harmonic 
expansion for r > 11 a, 

V harm (r) = U° + y (r- r 0 , h ) 2 - (8) 

We derive the mesoscopic parameters V^,fc^,ro,^ by fit¬ 
ting Vjiarm ) in a neighborhood of the minimum to the 
data V r (r) obtained from microscopic simulations. In 
figure 5 the resulting parameters and the two curves are 
shown. 

Moreover, we have compared V r (r) with the following 
expression for a finitely extensible non-linear elastic po- 
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Figure 4: Average distance between the two mesoscopic parti¬ 
cles versus the angles 6 \ and O 2 shown as a color and contour 
map for the two cases 0 = 0 (top) and 0 = 7r (bottom). It can 
be seen that even by rotating only one of the two particles, 
the average distance can be reduced considerably. However, 
the maximum reduction is observed when both particles are 
rotated. The variation of the torsion angle 0 does not change 
the main trend but has a moderate influence for intermediate 
rotations. 


tential (FENE) [42, 56], 


K f rl 


Vfene (r) = Vp -- 


-In 


1 - 


r - r 0 


(9) 


Figure 5: Effective pair potential V r (r) obtained from the 
microscopic simulation data and fit using a simple expression 
Vharm(r ) as in (8)- The fit was made in the interval [Her, 27 a]. 
In this way, the values for the mesoscopic model parameters 
ro,h,kh,Vh are determined. 



r/cr 


Figure 6: Effective pair potential V r (r) obtained from the mi¬ 
croscopic simulation data and fit using Vfene(v ) from (9) 
leading to the parameter values as listed in the plot. The 
fit was made in the interval [llcr, 52<r|* In this way, the val¬ 
ues for the mesoscopic model parameters r°, Vmax^Kf , Vp are 
determined. 


It takes into account that the chain cannot extend be¬ 
yond a maximal length, since Vfene (r) diverges when 
r -A ro + Vjnax. We see that (PVfene^)/ dr 2 = Kf for 
r = ro and therefore Kf represents the elastic constant 
in a harmonic expansion of this non-linear potential. As 
for the harmonic expression, we fit Vfene{ r) to our mi¬ 
croscopic data and thus derive the mesoscopic model pa¬ 
rameters Kf, ro, Vmax, and Vp as displayed in figure 
6. The agreement between the resulting expression for 
Vfene{ r) and V r (r) in the regime r > llcr is excellent, 
especially for the branch of the curve right to the min¬ 
imum. According to the fit, the maximum extension of 
the chain occurs for r = ro + r max — 75 a. In fact, since 
the radius of the mesoscopic particles is 5 a and each of 


the 60 beads making up the polymer chain has diameter 
cr, at r = 70 cr the polymer chain is completely stretched. 
A further elongation is of course possible due to the har¬ 
monic inter-bead interaction and this justifies the result 
of ^ 75a for the maximal extension. 

Last, we compare the elastic constants kh and Kf ob¬ 
tained from the harmonic and FENE approximation, re¬ 
spectively, as listed in figures 5 and 6. The resulting val¬ 
ues of 0.0 UksT/cr 2 and 0.015fc^T/a 2 are in good agree¬ 
ment with each other. 

To find a mesoscopic model expansion for the effective 
pair potential Ve 1 ( 6 1 ) [VJ? 2 (#2) has a very similar behav- 
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Figure 7: Effective pair potential Vq^Oi) calculated from the 
microscopic simulation data and fit using the phenomenologi¬ 
cal mesoscopic expression Vd(0\) from (10). Resulting values 
for the mesoscopic model parameters are listed in the plot. 
The fit was made on the interval [0, tt\. 

ior], we compare it with the phenomenological expression, 

Vb(0i) = Vg + £>[cos(0i) - cos(6»o)] 2 (10) 

introduced in (3) of [30]. As before, we can derive the 
mesoscopic parameters V^, D, and 0o by fitting Vd(0 i) 
to the microscopic data represented by Vq^Qi). The re¬ 
sulting parameters and the comparison between the two 
curves are shown in figure 7. Although (10) does not per¬ 
fectly reproduce the one-variable pair potential Vq 1 (6 i), 
it appears as a reasonable approximation in the neigh¬ 
borhood of the minimum energy. Moreover, the rather 
flat behavior of Ve 1 (6 1 ) for small 0\ is well represented 
by Vd(0i), which is at leading order proportional to 0\ A . 

Finally, we want to find a mesoscopic expression to re¬ 
produce the effective pair potential V^(0) obtained from 
the microscopic data due to relative torsional rotations 
between the two particles around the connecting vector 
r. We fit the microscopic data V^(0) around the mini¬ 
mum with the expression 

V T (<j>) = V T ° + r[cos(0o) - cos(0)]. (11) 

It leads to the parameters and fit depicted in figure 8 and 
is quadratic at lowest order in </>, V T {(j)) ~ V® +r</> 2 . The 
parameter 0o is redundant and can be absorbed into V^P, 
but we leave it for reasons of comparison to the following 
( 12 ). 

It is interesting to compare V^(0) with the expression 
taken from (4) of [30], 

V'M) = V'° T + t'[cos(<#,) - cos(<^)] 2 . (12) 

The discrepancy between the two curves shown in figure 
9 is obvious. The reason becomes clear when we expand 



Figure 8 : Effective pair potential V^(</>) calculated from the 
microscopic simulation data and fit using V T ((/)) from (11). 
Resulting values for the mesoscopic model parameters are 
listed in the plot. The fit was made in the interval [0,0.87t], 



Figure 9: Effective pair potential V^(0) calculated from the 
microscopic simulation data and fit using V' T {<j>) from ( 12 ). 
Resulting values for the mesoscopic model parameters are 
listed in the plot. The fit was made in the interval [ 0 , 0 . 87 r]. 

(12) to lowest order in 0, V f r {(j>) — V'® + r'0 4 . This 
expression is quartic in the torsion angle </>, leading to 
the comparatively flat behavior in the region around the 
minimum. The last two comparisons suggest that the 
analytical form of the mesoscopic model pair potential 
as a function of the azimuthal angle, which acts as a 
torsion, can be optimized by using a quadratic form as in 
(11) instead of the quartic one implied by (12) at leading 
order around the minimum. 

Finally, we performed additional microscopic simula¬ 
tions to estimate how the microscopic system parame¬ 
ters affect the mesoscopic model parameters. In par¬ 
ticular, the influence of the mesoscopic particle radius 
a and the number of beads N of the connecting poly¬ 
mer chain was investigated. As a trend, we found that 
the single-variable potentials tend to become stiffer for 
shorter chains and for bigger mesoscopic particles (see 

























































appendix A). Moreover, we observe that the fit of the 
functional forms in expressions (9)-(ll) with the simu¬ 
lation data further improves for increasing chain length 
and decreasing size of the mesoscopic particles [57]. 

So far, we have discussed mesoscopic analytical expres¬ 
sions to approximate the one-variable effective pair po¬ 
tentials V a {a). In the following section, using a numerical 
fitting procedure, we will determine harmonic coupling 
terms that take account of the correlations between the 
mesoscopic variables. In this way, we will further develop 
and improve our approximation of the probability density 
p c ( 7 ) from the microscopic simulations. 


VII. BUILDING A COUPLED EFFECTIVE 
PAIR POTENTIAL 

Our goal is to describe the effective interaction between 
the two mesoscopic particles, coarse-graining the role of 
the connecting polymer chain. The entropic nature of 
the polymeric interactions provides a direct route to av¬ 
erage out the microscopic degrees of freedom and thus 
build an effective scale-bridged model. A natural way 
to proceed would be to find an analytical approximation 
for p c ( 7 ) and thus derive an analytical expression for the 
mesoscopic pair potentials in the spirit of (3). 

As a first approximation we may describe p c as the 
simple product of the one-variable probability densities, 
n a Pa(&)- We obtain 

P'appi 7) = exp { —/3 [Vw ca (r/2a) + Vfene(j") 

+ v d ( 0 1 ) + v D (e 2 ) + v t (4>)]}, ( 13 ) 

when we use (7). Analytical approximations for the pair 
potentials V a (a) were derived in (9)-(11) in section VI 
together with the fitting parameters listed in figures 6- 
8. Care must be taken for the term V r {r), which has 
to be substituted by Vwca{^/ 2a) + Vfene (r) to take 
account of the steric repulsion between the mesoscopic 
particles, see (1). (13) correctly describes some aspects 
of the system behavior. For instance, the steep varia¬ 
tion due to the WCA potential between the mesoscopic 
particles is well represented in this approximation. How¬ 
ever, this description would lead to a distribution with 
vanishing correlations between the mesoscopic variables. 
That corresponds to assuming them independent of each 
other, which omits some important physical aspects, see 
the wrapping effect in section V. 

To make a step forward and take account of correla¬ 
tions we multiply correction terms to the previous fac¬ 
torized approximation in (13), obtaining the expression 


the one hand, we can simply fit the original probability 
density (e.g. minimize the squared difference between p c 
and Papp ). On the other hand, we can follow a moment¬ 
matching approach, looking for a p app that has a corre¬ 
lation matrix as close as possible to the original one. We 
performed a mixed strategy by fitting the expression in 
(14) to the simulation data p c ( 7 ), using as a criterion 
for best initialization of the fit an outcome that as close 
as possible matches the correlations directly calculated 
from the simulation data p c ( 7 ). This fit was performed 
by minimizing the squared difference between p c { 7 ) and 
Papp ( 7 ) 5 using the Nelder-Mead algorithm [58] provided 
by the SciPy library [59]. In this procedure, normaliza¬ 
tion of the approximated probability density as in (5) is 
enforced. The best set of parameters £ and S found is 


€ 


/23.54cr\ 

| 0.321tt 
0.321t r 5 


V 0 J 


(15) 


/ 4.810 _3 /(j 2 
0.19/(77r 
0.19/(777 
\—1.8 10 _3 /(77T 


0.19/(777 0.19/(777 

1.21/77 2 1.54/77 2 
1.54/77 2 1.21/77 2 
1.25/77 2 1.25/77 2 


-1.810“ 3 /(777\ 

1.25/77 2 I 
1.25/77 2 
-0.32/77 2 / 
(16) 


resulting in the correlation matrix 


8 


app 


1 

-0.306 
-0.306 
\ 0.045 


-0.306 -0.306 0.045 \ 
1 0.027 -0.071 1 

0.027 1 -0.071 

-0.071 -0.071 1 / 


(17) 


It would be unrealistic to try to exactly reproduce all 
properties of p c { 7 ) through an analytical approximation. 
Nevertheless, using the resulting expressions from (14)- 
(16), we can take account of the strong anticorrelation 
Q r o 12 between r and #i ; 2 , which is the strongest and 
most important one in the system; compare ( 6 ) and (17). 
For the other elements p r ^, qq 1 >2< £, and ge 1 e 2 we then ob¬ 
tain stronger deviations. However, those correlations are 
smaller than Q r o 12 a f least by a factor 4 and therefore 
carry a smaller amount of information about the overall 
system behavior. 

As a total result and in analogy to (3), we obtain from 
(14) the optimized analytical expression V app to model 
the effective interaction between the mesoscopic parti¬ 
cles: 


V app {if) = VwcA{r/2a) + Vfene{t) + Vb(0i) (18) 

+ V D (0 2 ) + V T {4>) + k B Tj2~ aP (a - £*)(/? - fr). 

a/3 


Popp(7 ) « p'a PP (l) X exp [- (7 - £) • S • (7 — £)], (14) 

where the elements of £, a 4-components vector, and S, 
a 4 x 4 symmetric matrix, are free parameters chosen to 
match the original data. There are at least two possible 
numerical approaches to find the best £ and S values. On 


The corresponding expressions and values for the fitting 
parameters are given by (1), (9)-(11), (15), and (16) as 
well as figures 6 - 8 . Thus, the effective pair potential is 
divided into two parts: one-variable and two-variable po¬ 
tentials. The former are the analytical single-variable 
pair potentials derived in section VI, see ( 1 ) and (9)-(11), 
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Figure 10: A simplified sketch of the mesoscopic configura¬ 
tion in which the magnetic moments are assigned to the par¬ 
ticles. The depicted orientations a\ and a 2 , which identify 
the anchoring points of the polymer chain, correspond to the 
0i = 62 — 0 and 0 = 0 configuration. The distance between 
the centers of the mesoscopic particles here corresponds to 
r = 26a. In this configuration, the magnetic moments mi 
and m 2 are introduced parallel to each other and pointing 
along the connecting vector r. 


together with the diagonal a = 0 terms in the double 
summation of (18). The latter are the off-diagonal a j- /3 
terms and take account, to lowest order, of the corre¬ 
lations between different mesoscopic degrees of freedom. 
Correlations between r and 9±^ are the dominating ones, 
leading to such physical effects as the wrapping effect 
introduced in section V. 


VIII. IMPACT OF MAGNETIC DIPOLE 
MOMENTS 

We will now consider how magnetic interactions be¬ 
tween the mesoscopic particles modify the physics of the 
system, in particular the probability densities and other 
mesoscopic quantities. For this purpose, we assign to 
each particle a permanent magnetic dipole moment 
(i = 1 , 2 ), in the present work not going into the de¬ 
tails of how these could be generated. We introduce 
these magnetic moments in the state of highest probabil¬ 
ity density. At vanishing magnetic moments, this is the 
state of minimal effective energy of the mesoscopic sys¬ 
tem, i.e., following (3), the one that has maximum p c {l) 
over all the configurations 7 . The maximum occurs for 
(r M , Of 1 , O 2 1 , 4 > M ) — (26cr,0.0,0.0,0.0), and is displayed 
in figure 10. When the mesoscopic particles are in this 
configuration, we assign to them the two magnetic mo¬ 
ments mi and m 2 that are parallel to each other and 
point along the connecting vector r, i.e. in their orienta¬ 
tion of minimal magnetic energy, as depicted in figure 10 . 
To first order, this configuration should leave the angular 
orientations unchanged when increasing the magnetic in¬ 
teraction, but see also the discussion in section V. This 
is just one of the possible orientations that the moments 
could assume. However, such a configuration could be 
achieved with a certain probability, for instance, when 


the sample is cross-linked in the presence of an external 
magnetic field that aligns the magnetic moments [60-63]. 
The dipole moments mi and m 2 are assumed to have 
equal magnitude m = |mi| = |m 2 | and are rigidly fixed 
with respect to each particle frame. We measure the mag¬ 
netic moments in multiples of mo = yj 47 rfc#T( 2 a) 3 //io, 
where /io is the vacuum magnetic permeability. Then, for 
each 7 , we can calculate the magnetic dipole interaction 
energy between the two moments 

t r _ Mo mi • m 2 r 2 - 3(mi • r)(m 2 • r) ^^ 

Vm\i) — A 5 (f^J 

47 r r b 

Pom 2 —2 cos^T 1 ) cos^™) + sin(^j n ) sin^™) cos( 0 m ) 

4tt r 3 

0 ™ indicates the angle between m* and r, while 0 m is the 
angle between the projections of mi and m 2 on a plane 
perpendicular to r. These quantities can be expressed in 
the 7 variables once the orientations of mi and m 2 with 
respect to the particle frames are set by the protocol we 
described above. 

For non-magnetic polymer chains, the potential acting 
on the mesoscopic level V ( 7 ) is separable into a sum of 
magnetic and non-magnetic interactions V{~y) = ^( 7 ) + 
Vm( 7 )- Consequently, magnetic effects do not modify the 
contributions of the polymer chain to V c ( 7 ), as is further 
explained in appendix B. 

This corresponds to a factorization of the probability 
densities p c (l) and p m ( 7 ), where the latter is defined via 
the Boltzmann factor 


Prn ( 7 ) = exp [-pv m ( 7 )] (20) 


in analogy to (3). Therefore the total probability density 
becomes 


pin) 


e -/3[Vc(7)+Kn(7)] 

Z(m) 


Pc(7)Pm(7) 

Z(m) 


( 21 ) 


where Z{m) = f p c ( / y)Pm( / y)d~/ is the partition function 
describing the system for non-vanishing magnetic mo¬ 
ments. We can calculate averages on the system with 
magnetic interactions by 

( 22 , 


The single-variable marginal probability density p a (a ), 
a = r, # 1 , 02 , 0 , is again defined as £>( 7 ) integrated over 
all the 7 variables except for a. This is the same pro¬ 
cedure as described in section IV but substituting p c ( 7 ) 
with p{ 7 ). As can be seen from figure 11 , p r (r) shows an 
increase of the probability to find the particles close to¬ 
gether when m is increased. A peak builds up at small r 
because the magnetic energy tends to V m (r —> 0) —>• — 00. 
Although the magnetic spheres attract each other, a col¬ 
lapse is prevented by the WCA-potential, which becomes 
effective at r < 11 a and for r —>> 0 diverges faster to +00 
than the magnetic energy to —00. For m = m 0 , the pres¬ 
ence of a double peak in p r (r) could be connected to a 
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Figure 11: Nonvanishing permanent magnetic dipole mo¬ 
ments of the mesoscopic particles and their impact on the 
system. We plot the marginal probability density p r for the 
distance r between the particles for different values of the 
reduced magnetic moment m/mo. Inset: average particle 
distance (r) as a function of the reduced magnetic moment 

m/mo. 



Figure 12: Marginal probability density pe 1 for the an¬ 
gle 0 1 (very similar for 62 ) for increasing reduced magnetic 
moment m/mo • Inset: standard deviation of 0±, ae 1 = 

61 — (#i)) 2 ) (very similar for O 2 ) as a function of m/mo • 



Figure 13: Marginal probability density p$ for the angle </> 
describing the relative torsion between the particles, for in¬ 
creasing reduced magnetic moment m/mo • Inset: standard 

deviation of 0, a$ — (((f> — (0)) 2 ) as a function of m/mo■ 

for any fixed 6 \ the minimum of Vm , and therefore the 
maximum of p m , is located at 62 = #1 and vice versa. 
As a consequence, the #1 = #2 = 0 configuration leads 
to the same magnetic interaction energy as 61 = 62 = tt. 
Concerning only magnetic interactions, both configura¬ 
tions show the same probability p m . Integrating out r, 
62 •> and <p fromp m ( 7 ), the shape of the resulting magnetic 
probability density for #1 is symmetric around 0 1 = 7 t /2 
with one peak at 6fi = 0 and one at 6fi = n. Therefore, 
coupling magnetic (p m ) and non-magnetic ( p c ) contribu¬ 
tions, we find that some probability shifts to higher val¬ 
ues of Ox due to the magnetic interactions. However, we 
find the standard deviation of pe 1 ( 0 i) to decrease with 
increasing m (see the inset of figure 12), meaning that 
the particles become less likely to rotate along the #1 (or 
likewise 62 ) direction. 

The same calculation for shows that the particles 

also become less likely to rotate around the connecting 
vector r. The probability density of </> ~ 0 rises and 
sharpens, as we can see in figure 13, and the standard 
deviation (see the inset of figure 13) decreases, confirming 
quantitatively the sharpening of p^(0). 



hardening transition of the kind described in [30]. Due 
to the mutual magnetic attraction between the parallel 
dipoles we expect the average distance ( r ) to decrease 
with increasing m, and indeed it does so, as can be seen 
from the inset in figure 11. 

The changes in the angular distributions for 6 \ and 62 
due to the magnetic interaction are illustrated in figure 
12. For the two angles 6fi and 62 the distributions po 1 ( 61 ) 
and po 2 (# 2 ) are similar, and the behavior for varying m is 
approximately the same. Therefore, we only display the 
results for po 1 ( 6 1 ). The magnetic moments tend to align 
in parallel along the connecting axis r, corresponding to 
their absolute energy minimum. Since we introduced the 
magnetic moments in that configuration for 6 \ =62 = 0, 


IX. THERMODYNAMIC PROPERTIES 

Finally, we provide the connection to the thermody¬ 
namics of our canonical system and demonstrate the in¬ 
fluence of the magnetic interactions. With the partition 
function Z(m) as in (21), we can calculate the overall free 
energy as F(m) = — &BTln[Z(ra)], the internal energy of 
the system as 

U(m) = (V)(m) = / F <T)Pc(7)P,n(7)dT 

Z{m) 

with V{^) defined before (20), and the entropy as 
S{m) = [U{m) - F(m)\/T. 
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Figure 14: Free energy F(m) = — ksThi[Z{m)\ as a func¬ 
tion of the reduced magnetic moment m/mo. The partition 
function Z(m) is calculated once using p c (~f) and once using 
Papp( 7 ). The result from the microscopic simulation is labeled 
with “T c ” and the one from the analytical approximation with 
“T ap p”. Note that for m = 0 the free energies are equal. 



Figure 15: Internal energy U(m) as a function of the reduced 
magnetic moment m/mo, calculated according to (23). The 
curve labeled with U U C ” shows the result from the microscopic 
data using (23), whereas in the one labeled with u U app ” Pc(t) 
has been replaced by p a pp(~i)- 


To test the validity of our effective potential scheme us¬ 
ing the coupling approximation described in section VII, 
we calculate the thermodynamic quantities using both 
probability densities pd'y) and p app ( 7 ) and compare the 
results. At vanishing magnetic moment, the free energy 
of the system is the same for both probability densities. 
This is expected, because it is a direct consequence of the 
normalization condition: Z(m = 0) = Z c = f p c ( 7 )d 7 = 
f Pa PP { 7)^7 = 1 by construction of p app ( 7 ). For non¬ 
vanishing magnetic moments, the partition functions ob¬ 
tained from the two probability densities (and thus the 
corresponding free energies) start to deviate from each 
other because the integral f p c {l)Pm{l)^l is, in gen¬ 
eral, different from f p apP (l)Pm(l)&l • With increasing 
magnetic interaction, the difference between p c ( 7 ) and 
Pa PP (l) becomes more important and, as shown in figure 
14, leads to a slightly increasing deviation in the free en¬ 
ergies F(m) calculated in both ways: at m ^ 1.224rao 
they already differ by ^ 5.6%. Analogously, we can com¬ 
pare the internal energies of the system, shown in figure 
15: at m = 0 the error due to the probability density 
approximation is ~ 0 . 12 % of the exact value, rising up 
to ~ 0.32% at m ~ 1.224rao. A similar deviation follows 
for the entropy, see figure 16, where, however, the error 
at m ~ 1.224rao increases to ~ 0.26% of the exact value. 
Overall, however, the relative errors remain small, which 
confirms the validity and viability of our coarse-grained 
effective model potential in (18). 

To summarize the effect of the magnetic dipoles on the 
system, we may conclude that the particles are driven to¬ 
wards each other. In other words, the average distance 
between them decreases (see figure 11 ) and the distri¬ 
butions for the particle separation and for the angular 
degrees of freedom sharpen (see figures 11, 12, and 13). 
This reflects a decrease in entropy, which becomes pos- 



m/rriQ 


Figure 16: Entropy S(m) = \U(m) — F(m)]/T as a function 
of the reduced magnetic moment m/mo calculated once using 
p c ( 7 ) from the microscopic simulation (labeled with “ < S' C ”), 
and once replacing p c {^() by p a pp(~i ) (labeled with u S apP ”). 


sible due to the gain in magnetic interaction energy. In¬ 
deed, the calculated entropy decreases with increasing 
magnetic interactions (see figure 16), which is achieved 
by the decreasing free energy and internal energy (see 
figures 14 and 15, respectively). 


X. CONCLUSIONS 

Most of the previous studies on ferrogels describe the 
polymer matrix as a continuous material [23, 32, 34, 35] 
or represent it by springs connecting the particles [30, 
36, 41-44], but only a few resolve explicitly the poly¬ 
meric chains [45, 46]. In particular, a link between such 
microscopic chain-resolved calculations and the expres¬ 
sions for the investigated mesoscopic model energies has 
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so far been missing. We have outlined in the present 
work a way to connect detailed microscopic simulations 
to a coarse-grained mesoscopic model. 

This manifests a step into the direction of bridging the 
scales in material modeling. Starting from microscopic 
simulations considering explicitly an individual polymer 
chain connecting two mesoscopic particles, we specified 
effective mesoscopic pair potentials by fitting analytical 
model expressions. In this way, we were able to opti¬ 
mize a coarse-grained mesoscopic model description on 
the basis of the input from the explicit microscopic simu¬ 
lation details. Furthermore, we have shown that correla¬ 
tions between the mesoscopic degrees of freedom must be 
taken into account to provide a complete picture of the 
physics of the system. Moreover, we have examined the 
effect of a magnetic interaction, finding a tightening of 
the system by reducing the average distance between the 
magnetic particles and by reducing the rotational fluc¬ 
tuations around the configuration of highest probability 
density. 
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Our system consisted of only two mesoscopic particles, 
for which we derived the corresponding effective pair in¬ 
teraction potential. However, using this pair potential in 
a first approximation, two- and three-dimensional struc¬ 
tures can be built up, similarly to elastic network struc¬ 
tures generated using pairwise harmonic spring interac¬ 
tions between the particles [36, 43]. Including as a first 
approach magnetic interactions between neighboring par¬ 
ticles only, the different angles between the magnetic mo¬ 
ment of a particle and the anchoring points of the poly¬ 
mer chain to its different neighbors must be taken into 
account. Yet, using the reduced picture of pairwise meso¬ 
scopic interactions, it should be possible to reproduce for 
example previously observed deformational behavior in 
external magnetic fields for two- and three-dimensional 
systems [45, 46]. 

The scope of our approach is a first attempt of scale¬ 
bridging in modeling ferrogels and magnetic elastomers. 
Naturally, the procedure can be improved in many differ¬ 
ent ways. For example, we would like to study a system 
with multiple chains connecting the magnetic particles, 
with anchoring points randomly distributed over the sur¬ 
faces of the particles. Such a development would elim¬ 
inate artificial symmetries in the model and be another 
step towards the description of real systems. Further¬ 
more, to study more realistic systems, interactions be¬ 
tween more than two mesoscopic particles would need to 
be considered. Likewise, also the effect of interactions 
between next-nearest neighbors connected by polymer 
chains could be included. The last two points imply a 
step beyond the reduced picture of considering only effec¬ 
tive pairwise interactions between the mesoscopic parti¬ 
cles. Finally, via subsequent procedures of scale-bridging 
from the mesoscopic to the macroscopic level [49], a con¬ 
nection between microscopic details and macroscopic ma¬ 
terial behavior may become attainable for magnetic gels. 


Figure 17: Effect of varying microscopic system parameters 
on the resulting mesoscopic model parameters. In the micro¬ 
scopic simulations, we varied the mesoscopic particle radius 
a and the number of beads N forming the connecting poly¬ 
mer chain. Fits of the single-variable mesoscopic model po¬ 
tentials (9)—(11) to the microscopic simulation data lead to 
the presented mesoscopic parameter values, (a-d, O) Model 
parameters for the FENE potential (9); (e,f, A) model pa¬ 
rameters for the angular 0 1 - and ^-potential (10); (g,h, □) 
model parameters for the torsional ^-potential (11). Values 
for 6> 0 and (po vanished in all cases and are not shown. The 
data points for a/5cr — 1 and iV/60 as 1 correspond to the 
results presented in figures 6-8. 
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Appendix A: Dependence on microscopic system 
parameters 

In sections VI and VII we derived an analytical ap¬ 
proximation for the effective mesoscopic model potential 
describing our system. Here, we test how the mesoscopic 
model parameters depend on the microscopic system pa¬ 
rameters. In particular, we performed additional micro¬ 
scopic simulations with varied radius a of the mesoscopic 
particles and varied number N of beads forming the con¬ 
necting polymer chain. 

In each case, we repeated the analysis of sections IV, 
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Figure 18: Elements of the correlation matrix g for differ¬ 
ent values of a and N. The data points for a/ha = 1 and 
iV/60 = 1 correspond to the results presented in (6). 


VI, and VII by fitting the mesoscopic single-variable 
model potentials (9)—(11) to the microscopic simulation 
data. This reveals the trends in the dependences of the 
coefficients in the mesoscopic model potentials on the 
parameters a and V, as depicted in figure 17. From the 
trend of the resulting model parameters Kf,D, and r in 
figure 17 (c,e,g) we conclude that the potentials tend to 
become stiffer as the chain becomes shorter or - at least 
for the rotational degrees of freedom - the mesoscopic 
particles become larger (see [57] for corresponding fitting 


curves). 

Moreover, we have calculated the trend in the corre¬ 
lations between the mesoscopic degrees of freedom for 
varying values of a and V, see figure 18. We find that the 
magnitude of the correlations decreases with increasing 
N or decreasing a. Thus, quite intuitively, the coupling 
between the 7 variables tends to decrease with longer 
chains or smaller mesoscopic particle sizes. 


Appendix B: Separability of the Hamiltonian and 
Consequences for Coarse-Graining 


In section VIII we have included the magnetic inter¬ 
actions analytically on the mesoscopic level. They had 
not been included in the microscopic simulations. This 
procedure is possible due to a separability of the mag¬ 
netic and non-magnetic effects which as a consequence 
implies a factorization of the corresponding probability. 
Thus the contribution of the polymer chain needs to be 
simulated explicitly. The contribution of any interaction 
acting solely on the mesoscopic particles can be exactly 
taken into account separately afterwards. 

This argument relies on the separability of the Hamil¬ 
tonian into a sum of mesoscopic and microscopic parts, 
as well as on the fact that magnetic interactions affect the 
mesoscopic part only. We write down the Hamiltonian of 
the system as 


= Tmesoh) + V^’ B ( 7 ) + V$* A { 7 ) 


(Bl) 

N 

N 


+ TmicroA) + £ Ww’Sa( 7, T) + T) 

+ E V wcaA) 

(B2) 

n= 1 

n<n' 


N—l 



+ v ^\ 1 ,t) + v^ b ( 1 ,t)+Y. v H {n+1 \ r )- 


(B3) 


n=1 


Similarly to the main text, let us here indicate with 7 
the degrees of freedom (now velocities included) of the 
mesoscopic particles and with T the ones of the micro¬ 
scopic particles that build up the chain. For brevity we 
here label the mesoscopic particles by A and B and the 
microscopic ones by the discrete indices rz, ml = 1,..., N. 
Moreover, we denote by ( 7 , T) the corresponding in¬ 
teraction between particles p and g, where the explicit 
expressions of Vwca, Vh, and V m are as described in 
(1), (2), and (19). Tmeso(l) and Tmi C ro<X) indicate, re¬ 
spectively, the kinetic energies of the mesoscopic and mi¬ 
croscopic particles. Therefore, 7 /( 7 , T) is separable into 


where 'H rneso ( 7 ) contains the terms in line (Bl) and 
Bmicro( 7 ) is composed of the terms in lines (B2) and 
(B3). 

Since we work in the canonical ensemble, the physics 
of the system derives from the partition function 


Z{m) = I e- /3w( ' |, ’ r) d 7 dr (B5) 

q — PT-Lmesoilf) e — (3'Hmicro('lf,r) 


7/(7, r) = l~Lmeso (7) + Umicroil, r), (B 4 ) with f3 = l/ksT. Coarse-graining means to integrate out 
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the microscopic degrees of freedom T, so we rearrange 
Z(m) = J e-P^rnesoh) ^ PUmicro foT) qp ^ (B6) 

— / p — f3'Hmeso(lf) p' (^VW 

/ c "micro \ I I 

= J e-P v ™ B ^e~P{ Tmeso ^ +v wcA^\z m i cro {~i) d7. 

The connection to the probability density p c { 7 ) is given 
by 


Pc (7) = 


^[lmeso(7) + hrCAd)] 

Z(m = 0) 


(B7) 


Since the magnetic interactions only affect the meso¬ 
scopic particles, see line (Bl), it is solely contained 
in Hmesoi y)- Therefore, the microscopic Hamiltonian 
'Hmicro (y, T) for a specific fixed configuration 7 of the 


mesoscopic particles does not explicitly depend on the 
magnetic moments. The magnetic effects do not modify 
the physics of the polymer chain for a given configuration 
of the mesoscopic particles. The underlying physical rea¬ 
son is that the polymer chain does not contain magnetic 
components. 

As we can see, the magnetic interactions only appear 
on the mesoscopic level of the final partition function 
Z. The magnetic interactions are simply included by 
multiplying on the mesoscopic level with the additional 
probability factor ex.p[—/3V r ^ ,B ('j)\ The remaining part of 
the integrand that contains the information drawn from 
the MD simulations, is not affected. 

Explicitly introducing the magnetic moments in the 
microscopic simulation would therefore not affect the fi¬ 
nal outcome. Thus, it is sufficient to determine the effects 
of the polymer chain through MD simulations and add 
the magnetic interactions analytically in a subsequent 
step. 
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